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We consider astrophysical maser radiation that is created in the presence of 
mildly supersonic, magnetohydrodynamic (MHD) turbulence. The focus is on 
the OH masers for which the magnetic field is strong enough that the separa- 
>■ I tions of the Zeeman components are greater than the spectral linebreadths. A 

^ ' longstanding puzzle has been the absence of the Zeeman vr components and the 

^ ' high circular polarization in the observed spectra of these masers. We first argue 

that the elongation of eddies along the field that has recently been recognized 



o 



in MHD turbulence will enhance the optical depth parallel to the magnetic field 
in comparison with that perpendicular to the magnetic field. We then simu- 
Qh' late maser emission with a numerical model of MHD turbulence to demonstrate 

O ' quantitatively how the intensities of the linearly polarized vr components are sup- 

^ • pressed and the intensities of the nearly circularly polarized a components are 

■ enhanced. This effect is also generic in the sense that most spectral lines in MHD 
turbulence with Mach number M ~ 1 should have larger optical depth parallel to 

^ ' the magnetic field than perpendicular. The effect is reduced considerably when 

■ M < 1. The simulations also demonstrate that the velocity and magnetic field 
variations due to the turbulence can (but do not necessarily) cause one of the 
a components to be much more intense than the other, as is often observed for 
mainline OH masers. 

Subject headings: masers: OH — polarization — magnetic fields — MHD — turbulence 



^Department of Physics. 

^Permanent address:Institute of Astronomy of the RAS, 48, Pyatnitskaya str., 109017 Moscow, Russia. 
■^Department of Astronomy. 



- 2 - 



1. Introduction 

Magnetohydrodynamic (MHD) turbulence is highly anisotropic (Goldreich & Sridhar 
1995, hereafter GS; see also Montgomery & Turner 1981; Higdon 1984; for numerical con- 
firmation see Cho & Vishniac 2000 and Maron & Goldreich 2001). According to the GS 
theory, turbulent eddies are elongated parallel to field lines, and this anisotropy increases 
as one progresses to smaller and smaller scales. In traversing a turbulent MHD gas, spec- 
tral line radiation will then typically sample regions that are extended parallel to the field. 
That is, the optical depth for spectral lines in MHD turbulence will be larger parallel to the 
field than perpendicular to the field. The GS theory predicts anisotropic scattering of radio 
waves, consistent with observations (e.g. Narayan, Anantharamaiah, & Cornwell 1989). It 
also predicts an angle-averaged spectrum that is consistent with measurements of radio wave 
fluctuations in the local interstellar medium (Armstrong, Rickett, & Spangler 1995). 

Radiative transitions of atoms and molecules in the presence of a magnetic field are 
split into a and vr components according to the change in the quantum number for angular 
momentum along the magnetic field Am = ±1 or Am = 0, respectively. The absorption and 
emission of radiation in these transitions by an isolated atom or molecule depends upon the 
angle 9 between the magnetic field and the direction of propagation of the radiation, and is 
proportional to (1 -|- cos^ 9) for a transitions and to 2 sin^ 9 for tt transitions (e.g., Condon & 
Shortley 1970, pg. 99). Radiation emitted in a transitions is elliptically polarized — ranging 
from entirely circular at ^ = to entirely linear at ^ = 90°, whereas that due to vr transitions 
is 100% linearly polarized at all angles. Hence, if the magnitudes of the optical depths 
in directions close to that of the magnetic field are larger than in other directions, the a 
components and the circularly polarized radiation will tend to be enhanced in comparison 
with the TT components and the linearly polarized radiation. For masers, the effect of such 
differences can be exponentially amplified. 

Among interstellar spectral lines, only the OH masers have Zeeman splittings that are 
greater than the spectral line breadths and they do, in fact, exhibit the characteristics de- 
scribed above. The OH masers at 1665 and 1667 MHz ("mainline" masers) occur widely in 
regions of star formation and have been especially well studied (e.g., Garcia-Barreto et al. 
1988; Argon, Reid, & Menten 2000). Except possibly for a recent observation (Hutawarakorn, 
Cohen, & Brebner 2002), the tt components have never been detected in this maser radia- 
tion. Circular polarization dominates and in some regions the radiation is essentially 100 
% circularly polarized. Nevertheless, a small fraction of the radiation often is linearly po- 
larized. These spectral characteristics have remained a puzzle (e.g., Reid 2002) since the 
initial discoveries of masers in astronomy. Though the Zeeman splitting ordinarily is much 
smaller for other masing transitions of OH, the tt components seem to be absent and circular 
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polarization also seems to dominate for these (Baudry & Diamond 1998 [13 GHz] ; see also 
Baudry et al. 1997 [6 GHz]). Gray and Field (1995) have emphasized that the a and tt 
radiation of the OH masers will tend to be beamed parallel and perpendicular, respectively, 
to the magnetic field. 

A secondary polarization characteristic of the OH masers is that, commonly, only one 
of the two cr components of the 1665/1667 MHz transitions is detected from a particular 
maser location in observations that are performed at the highest possible angular resolution, 
though comparable numbers of masers with single right or single left circular polarization 
may be observed from the entire cluster of masers. For other masing transitions of OH, 
both a components more often are detected at the same location (e.g., Moran et al. 1978; 
Caswell & Vaile 1995). The tendency for one of the a components of the 1665/1667 MHz 
transitions to be much stronger than the other has often been attributed to shifts in the 
resonant frequency caused by gradients in the Doppler velocity and in the magnetic field that 
partially cancel for one of the a components, but not for the other (Cook 1966). Actually, a 
velocity gradient alone is sufficient to create a single, sharp spectral line with a single sense 
of circular polarization if the change in velocity through the maser is as large as the Zeeman 
splitting and the maser is at least partially saturated (Deguchi & Watson 1986; Nedoluha 
& Watson 1990). In either case, the observation that single components with both senses of 
circular polarization occur in close proximity is more suggestive of the randomness associated 
with turbulence than to large scale gradients, which would tend to produce only one sense 
of the circular polarization for masers that are located close to one another. 

The calculations of this Paper explore in more detail how well the anisotropy of the 
medium caused by MHD turbulence can serve as the basis for an understanding of the po- 
larization characteristics of the OH masers, and conversely, whether these observed charac- 
teristics may be evidence of the nature of MHD turbulence in the interstellar gas. Standard 
methods are employed for computing a number of simulations of turbulent, compressible 
MHD to serve as representative "cubes" of velocities, magnetic fields and densities of the 
masing medium. These cubes are used to represent the regions of interstellar clouds in which 
clusters of OH masers are found. The coupled radiative transfer equations for the Stokes 
intensities are then solved as a function of Doppler velocity by numerical integration for a 
large number of uniformly spaced rays of maser radiation that pass through the cube at 
various angles. A comparison of the optical depths computed without the effect of maser 
saturation for the rays parallel and perpendicular to the direction of the average magnetic 
field provides a robust indication that the anisotropy of MHD turbulence has a significant 
effect. In examining the computed spectra, we focus on the brightest rays. These cover only 
a small fraction of the entire surface, just as the observed masers cover only a small fraction 
of the surface of the cloud in which they occur. We begin with general considerations for 
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both the MHD turbulence and the maser polarization (§2) before giving a more detailed 
description of the basic equations for the calculations (§3). The results of the computations 
are presented in §4. In §5, we discuss further the requirements for the MHD anisotropics 
to be applicable in the environments of the OH masers. We also reason in §5 that Faraday 
rotation, which has been mentioned as a possible cause for the absence of the tt components 
of OH masers, is an unlikely explanation. 



2. General Considerations 

2.1. MHD Turbulence 

The statistical properties of MHD turbulence can be partially characterized by the 
structure function 

vl), = {(v,(r)-v,(v + X)r). (1) 

We will use the shorthand to indicate the structure function for the velocity component 
parallel to the magnetic field, and v± to indicate the structure function for either of the two 
perpendicular components. The quantity vf^ means {v^^ + 2v\) averaged over all directions 
for the separation vector A. For Kolmogorov-like turbulence, vx ~ vl{\/ L)^/^ , where vl is 
the three dimensional velocity dispersion on the outer scale L. 

At any scale A it is useful to consider two characteristic dimensionless parameters: 
the sonic Mach number Mx = vx/cg and the Alfven Mach number Ma,x = vx/va (we will 
assume that p and va do not scale with A) . When Ml = M <S 1 turbulence is approximately 

incompressible; when Ma <^ 1 it is strongly constrained by the presence of the mean magnetic 
field, which it cannot bend. Because vx is typically an increasing function of A, Mx and 
Ma,x ai'e both small at small A. MHD turbulence is therefore incompressible and mean field 
dominated on small scales. This was the case considered by GS in their theory of strong 
MHD turbulence. 

The OH maser lines, however, will be little infiuenced by turbulence unless the Doppler 
shift associated with bulk motion of the gas is at least comparable to the thermal linewidth 
of OH. This requires Ml > 1, which implies that the turbulence is compressible. There is 
currently no good theory for compressible MHD turbulence. Simulations (Cho & Lazarian 
2003; Vestuto, Ostriker, & Stone 2003) indicate, however, that the scahngs introduced by 
GS may be valid up to Ml ~ a few. Encouraged by this, we will employ GS's scalings in 
the estimates that follow. 

We want to estimate the line-center optical depth tq in an unsaturated maser line. For 
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simplicity we will make an estimate assuming that the optical depth at Doppler velocity v 
has the form 

T^^ — ! dse-^''-''*"''''^'^^'/'''^ (2) 

Vth J 

where cr is a constant, Vturb{s) is the component of the turbulent velocity along the line of 
sight, Vth is the thermal velocity of OH, and the integral is taken along the line of sight. Most 
of the optical depth will come from a region of size Ifh where Vturb{s) changes by of order 
Vth- If kh is smaller than the size of the masing cloud (i.e., vl > Vth), then tq ~ {cr/vthjlth- 

Otherwise we simply have tq ~ {a/vtii)L. We will assume that we can estimate Ith using the 
structure function (1). Thus to find the line center optical depth we need to find kh such 
that vi^^ - Vth- 

Perpendicular to the mean magnetic field, the GS theory gives 



as in the Kolmogorov theory. Then 

V^Vth 



Hh,± 



Vl 



(4) 



where we have normalized the perpendicular velocities by assuming isotropy at the outer 
scale: 

-II - <L - f ■ (5) 

A more general treatment is possible, but this simplified analysis will serve to make our 
point. 

To find lth,\\ we must invoke a weak extension of the GS theory, because GS theory 
considers incompressible Alfvenic fiuctuations which are necessarily perpendicular to the 
mean field. But models of weakly compressible turbulence (Lithwick & Goldreich 2001), 
numerical models of compressible turbulence (Cho & Lazarian 2003; Vestuto, Ostriker, & 
Stone 2003), and of the cascade of slow (pseudo-Alfven) waves in incompressible turbulence 
(Cho, Lazarian, & Vishniac 2002) suggest that v\\ should have the same spectrum as v±. 

Correlations along the field can be determined from GS's relation Ay ~ A^''^, which 
defines a surface of constant correlation amphtude. Then 



1/2 



Mh)^^(^) , (6) 
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and 




2 



(7) 



Again, we have assumed isotropy at the outer scale. 



Then the ratio of parallel to perpendicular optical depths is 

kh,\\ _ Vl 



(8) 



kh,± V^Vth ' 



provided that vl > VthV^ (if the turbulence is subthermal throughout the emitting region 
then the parallel and perpendicular optical depths differ only slightly). Parallel optical depths 
should therefore typically be larger than perpendicular optical depths for mildly supersonic 
turbulence. 

When will this effect be absent? The parallel optical depth can exceed the perpendicular 
optical depths only if vi/vth ^ \/3- Notice that vi/vth = ML^/rna^j2ii pa 2.1Ml, where n 
is the mean molecular weight of the gas and mon is the mass of the OH molecule, so only 
Ml ~ 1 is required for an interesting ratio of optical depths parallel and perpendicular to 
the field, and conversely the effect should vanish for subsonic turbulence. 

The key is, of course, the anisotropy of the turbulence. We expect that MHD turbulence 
is nearly isotropic unless Ma,i^,^ < 1. That is, the field needs to be strong enough to place a 
dynamical constraint on the masing eddies. Since 



The estimates in this section suggest that the optical depth to OH maser lines should 
be larger parallel to the mean magnetic field than perpendicular to it, provided that the 
turbulence is mildly supersonic and (3 <1. Our estimates are encouraging but not conclusive 
because they rely on an extension of the GS theory. In addition, strong maser spots are rare 
events whose frequency may not be well described by the two-point statistics considered 
here. A test of these estimates against a simulation of compressible 3D MHD turbulence is 
required. 




(9) 



the effect is absent unless P <1. 



2.2. 



Polarized Maser Radiation 



The transfer equations for polarized radiation (see §3) can readily be solved for a uni- 
form medium in the regime where the maser is unsaturated and the Zeeman splitting of the 
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components is much greater than the spectral hnebreadth so that overlap can be ignored 
(Goldreich, Keeley, & Kwan 1973, hereafter GKK, pg. 124). For an optical depth of mag- 
nitude To, the intensities I± of the two a components and Iq of the tt component are given 

by 

/± ~ ( V2) exp[To(l + cos' e)/2] (10) 

and 

/o~(/e/2)exp[Tosin2^] (11) 

for propagation at an angle 9 from the magnetic field when the amplification is large as is 
the case for astrophysical masers. Here, Ic is the unpolarized continuum radiation that is 
incident on the far side of the maser and serves as the "seed" radiation. In terms of the 
Stokes parameters V and Q that measure the circular and linear polarization, respectively, 
the fractional polarizations are 

\4 _ ±2cose 

I± 1 + C0S2^' ^ ' 

Q± ^ sin' 9 

7± l + cos^e' ^ ' 

0, (14) 



h 
and 

^ = -1, (15) 

-'0 

where the direction of the linear polarization is parallel (perpendicular) to the direction of 
the magnetic field projected on the sky when Q is negative (positive). A reasonable estimate 
for the mainhne OH masers is Tq ~ 20. Prom equations (10) and (11), the intensity I± 
of the a components can then be seen to peak strongly at small angles 9 where circular 
polarization dominates. The intensity /q of the linearly polarized tt components is similarly 
peaked in directions that are perpendicular to the magnetic field. We would not then expect 
to observe the a and tt components together from the same region if the magnetic field 
there is uniform and the masers are unsaturated. However, statistically we are equally likely 
to observe regions where the magnetic field makes large angles with the line-of-sight and 
regions where it makes small angles. Hence, in conflict with the observations, comparable 
numbers of cr and tt components would be detected when a number of independent regions 
are observed. The situation is even worse for understanding the absence of tt components 
in the limit of highly saturated masing. All three intensities {I± and Iq) become equal and 
independent of angle 9 with fractional polarizations that are the same as given by equations 
(12)-(15) for unsaturated masers (GKK, pg. 122). 
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On the other hand, if the medium is anisotropic and the optical depth for unsaturated 
masing is greater in directions close to the magnetic field than perpendicular to the field 
lines, then I± near ^ = can be much greater than Jq near 9 = 90° (or at any other angle) 
even for modest fractional differences in the optical depths according to equations (10) and 
(11). It is then possible that the a components near ^ = are strong enough to be detected 
while the tt components near 9 = 90° or at any other angle are too weak to be detected. As 
described above, the anisotropy of MHD turbulence can cause the optical depths along the 
field hues to be greater than those perpendicular to the field lines. 

Trapped infrared radiation and perhaps other processes can randomize the thermal 
velocities and the populations of the magnetic substates of the masing molecules at a rate 7. 
This rate can be greater than the rate F at which an excited masing state decays as a result 
of all processes other than the stimulated emission by the maser radiation which proceeds 
at a rate it! (Goldreich, Keeley, & Kwan 1973; Anderson & Watson 1993). For simphcity, we 
take the rate 7 to be the same for randomizing both the velocities and substate populations. 
A regime 7 > i? > F then exists in which the maser is saturated in the sense that the 
maser radiation influences the molecular populations and the amplification is approximately 
linear (as opposed to exponential) with the length of the maser in a uniform medium, but still 
behaves like an unsaturated maser in retaining angular distributions for a and tt components 
of the form given by equations (10) and (11). The spectral line narrowing of unsaturated 
maser amplification to produce subthermal linebreadths also is preserved in this regime. Just 
as for completely unsaturated masers, the velocity distribution of the masing molecules is 
Maxwellian and the populations of an energy level are equal in this regime. 



Our analysis is based on a numerical model of decaying MHD turbulence. In the model, 
we integrate the equations of compressible, ideal MHD: 



3. 



Basic Methods 



3.1. 



Numerical Model 



dp 

at 



+ V • (pv) = 0, 



(16) 




dv 



Vp (B • V)B 



(17) 



dB 

'dt 



V X (v X B), 



(18) 
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and as usual p = density, v = velocity, p = pressure, and B = magnetic field. We use an 
isothermal equation of state p = cj.p, with = constant. This equation of state is valid in 
regions where the cooling time is much shorter than a dynamical time. No forcing term is 
employed, and the model is nonself-gravitating. The scheme contains no explicit resistivity 
or viscosity, although a nonlinear artificial viscosity is used to capture shocks. 

The basic equations are integrated in a cubic, periodic domain of linear size L. The 
initial state consists of a uniform medium, p = Po = constant, B = B^o^, with a superposed 
velocity perturbation 6v. The 5v is a Gaussian, random field with the following properties: 
(1) the spectrum (5v^) ~ k~^^^^ for kmax > k > A;mm(this is a Kolmogorov spectrum); (2) 
kmax = 32(27r/L); (3) kmin = 2(27r/L), (4) V-5v = 0; (5) the modes with k = kmin had their 
amplitudes set equal to the expectation value. 

The power-law spectrum was adopted to hasten the relaxation of the turbulent state. 
A Kolmogorov slope was chosen (rather than the —4 expected of strongly compressible 
turbulence) because our models have a Mach number of a few and are therefore only weakly 
compressible. Experience suggests that the evolution is insensitive to the precise choice of 
index as long as most of the power is concentrated in the long wavelength spatial modes. 

A fixed (rather than resolution-dependent) kmax was introduced to permit evolutions 
of identical reahzations of the power spectrum at different resolutions. The kmin was set 
at a wavenumber corresponding to a wavelength L/2 rather than L in order to increase 
the number of modes that contribute significantly to the power. This reduced run-to-run 
variations in the evolution, and was also the motivation for fixing the amplitude of modes 
with k — kmin- The initial velocity field was made incompressible to facilitate comparison 
with earlier runs (e.g. Stone, Ostriker, & Gammie 1998) which used a similar constraint in 
order to minimize turbulent dissipation. 

The initial conditions are scaled so that L — Cs — po — 1. These quantities can then be 
rescaled to match physical conditions in a masing region. 

Our numerical method is based on the ZEUS scheme (Stone & Norman 1992). ZEUS 
is an operator-split, finite difference scheme on a staggered mesh. Constrained transport 
(Evans & Hawley 1988) is used to guarantee that V ■ B = to machine precision. Our 
version of the ZEUS code was parallelized by one of us (JCM) and run on NCSA's Platinum 
cluster. 

Our numerical models have 3 parameters: the rms Mach number of the initial conditions, 
the mean magnetic field strength (which can be parameterized in terms of the initial f3~^ = 
v\l(?g = BIq/IAttpocD), and the numerical resolution. We have studied runs at resolution 
66^, 128^ (our default resolution) and 256^. Our results are insensitive to resolution. The 
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initial rms Mach number in all of our computations is four. Three values for the initial i'a/cs 
are considered — 1, 3, and 10. For each of these wa/cs, the MHD equations are integrated 
for three independent choices for the ensemble of random numbers that specifies the initial 
amplitudes of the various Fourier components ^v^. We thus obtain three sequences in time 
of MHD cubes to represent the interstellar gas for each va/cs- 



3.2. Maser Radiative Transfer for the Stokes Pcirameters 



The mascr radiation is calculated with the simplification that each ray is treated as 
an independent linear maser — a standard simplification for astrophysical mascrs. Since the 
focus in calculating the intensities is on only the small fraction of the rays through the 
computational cube that arc the brightest and represent the observed maser features, the 
assumption that they can be treated as independent is reasonable since it is less likely 
that they will intersect the paths of other strong rays that are propagating in different 
directions. In any case, interaction between the rays can only have an effect when the maser 
is saturated and the effect would be to strengthen the stronger maser beam at the expense 
of the weaker. Since the cr radiation tends to be the stronger in the scenario here, the 
effect would be to further enhance the a radiation relative to the vr radiation, and hence to 
strengthen the conclusions of this investigation that MHD anisotropy tends to suppress the 
TT component. In addition, the present calculations are limited to masing under conditions 
of modest saturation. 

The calculations are performed for an angular momentum J = 1 ^ transition and, as 
is commonly done, the results are assumed to be indicative for molecular states with higher 
angular momenta. For the OH masers, it is clear that the Zeeman frequency (sometimes 
"gff2" ) is much larger than the rates for stimulated emission and other competing processes. 
In this regime, the radiative transfer equations for the Stokes parameters I,Q,U, and V are 
of the same form as for ordinary spectral line radiation and for the purposes here can be 
expressed as 

/ / \ 

Q 
u 

\v J 

from Watson & Wyld (2001) (which rephrases GKK for the relevant regime) where 

A={1 + cos^ 0)n{f+ + /_) + 2n/o sin^ 6, 
S = sin2e(/+ + /_-2/o)n, 





( ^ ) 
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B 
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U 











(19) 



(20) 
(21) 
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and 



C = 2ncose(/+-/_). 



(22) 



Here, the fs are the MaxweUian distributions for the component of the molecular velocity 
along the line-of-sight and n is the difference between the normalized populations of the 
magnetic substates of the upper J = 1 energy level and the normalized population of the 
J — lower energy level. The f± and /o are evaluated at the velocities v± and Vq obtained 
from (1 — v±/c)u = Ufi± and (1 — Vo/c)uj = uir where oj is the angular frequency 

of the intensities in equation (13), uor is the resonance frequency of the masing transition 
for a molecule at rest and ±(?0/2 are the splittings of the m = ±1 substates due to the 
Zeeman effect. Only a single population difference enters in equations (20)-(22) because the 
calculations will be restricted to the regime of saturation (7 > Ft) where the populations of 
the magnetic substates of an energy level are equal due to the assumed rapid cross/velocity 
relaxation. The population differences n are found by solving the rate equations in this 
regime using the standard idealization of "phenomenological" pump A and loss F rates 



Here and Rq are the rates for stimulated emission for the m = ±1 and substates 

(see Watson & Wyld 2001 for detailed expressions for and Rq), and AA is the 

difference between the pump rates per substate into the upper and lower energy levels. Also, 
note that in the above definitions AA does not involve the MaxweUian velocity distributions 
which are included as the / 's. As a further simplification, only the values of 
and i?o at the line centers of the a and tt components in the local rest frames moving 
with the turbulent velocities are used in equation (23). The completely accurate procedure 
would involve averaging over the absorption profile. The inaccuracy introduced by this 
simplification is negligible for this investigation since the emergent intensities near line center 
are adequate for the purposes here, and since the degree of saturation will be modest (and 
we have confirmed that the differences between the results with and without saturation are 
small) . 

External continuum radiation is assumed to provide the seed radiation for the maser 
in equation (19). Whether external continuum or spontaneous emission provides the seed 
radiation is unimportant for the emergent radiation. The distance s is measured along 
the straight-line path of the ray and the Stokes intensities in equation (13) are expressed 
in dimensionless form. They are actual intensities divided by a characteristic "saturation 
intensity" Ig = Shw^T/lSTTC^AEAQ), where Ae is the Einstein A-value for the transition 
and AQ is the solid angle into which the maser radiation is beamed. When 1 = 1 the rate 
for stimulated emission is approximately equal to the decay rate T for the molecular states. 



n = 



AA/r 



(23) 



l + 4(i?+ + i?_ + i?o)/3r 
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The uncertain beaming angle of the maser radiation is thus incorporated into the definition 
of Ig as is customary for relating the intensity and the mean intensity "J" for the linear 
maser (e.g., Watson & Wyld 2003). 

Equation (19) is expressed in the coordinate system in which the axis of quantization 
for the magnetic substates is along the magnetic field. That is, the "z-axis" is along the 
magnetic field. Stokes Q and U also depend on the orientation of the coordinate system and 
hence upon the direction of B. However, in a turbulent medium the direction of B changes 
along the path of a ray of maser radiation. It is thus necessary to transform Q and U with a 
rotation of the coordinate system at each step of the integration along a ray to express Q and 
U in the coordinate system in which equation (13) is applicable. If the magnetic field rotates 
by an angle (f) when projected onto the plane perpendicular the direction of propagation, 
Q and U in the coordinate system after the rotation are related to the same quantities Q' 
and [/', but expressed in the system before the rotation, by Q = Q' cos 20 + U' sin 2(f) and 
U = -Q'sin2(/) + C/'cos2(/) (e.g., Chandrasekhar 1960, pg. 34). 

When turbulent velocities along the line-of-sight Vfurb ^-re incorporated explicitly and 
gO, is expressed in terms of the magnetic field and the magnetic moment /i, the /'s become 
f± = Fex.p{-{v - vturb ± I^By/vth) and /o = Fexp(-(f - vturbf/vl^). In these, v and 
Vth are, respectively, the Doppler velocity at which the Stokes intensities in equation (19) 
are specified and the thermal velocity Vth — {2kT /moHY^'^ of the OH molecules (F is the 
normalization constant for a Maxwellian). 

To solve equation (19), it is necessary to specify (1) the radiation that is incident on 
the far side of the masing cube to serve as the seed radiation, (2) the Zeeman splitting (or 
more correctly /x), and (3) the pumping rate AA. 

The incident continuum radiation is assumed to be unpolarized so that the initial values 
for the Stokes intensities are Q — U = V — Q and / = Ic- For an assumed brightness 
temperature of 30 K for the external continuum seed radiation, decay rate V = 0.03 s~^, and 
a beaming angle Afl — 5 x 10~^ ster, Ic — 10~^° when expressed in units of Is- 

All of the computations in the Figures are obtained with this 1^ though we have verfied 
that the results are insensitive to changes in Ic by a few powers of ten. For the Zeeman 
splitting, we adopt {^Bavg/vth) = 9 for the computations in the Figures. This corresponds 
to a separation of 5.4 km s^^ between the two a components when the gas temperature is 100 
K. We have also performed computations with {^Bavg/vth) = 4.5 and have verified that the 
essential behavior of the spectra is unchanged within this range of Zeeman splittings which 
is representative of what is observed for the OH masers with distinct Zeeman components. 
We choose to specify /iBavg/vth instead of simply /i so that in all of our computations, the 
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Zeeman components will be well separated as in the observations and the splitting will be 
essentially the same for purposes of comparison. 

We consider two ideahzations for AA: A A is equal to the density multiplied by a 
constant, and AA is simply constant throughout the computational cube. All of the results 
that we present in the Figures are obtained with AA proportional to density. The results 

are generally similar when AA = constant is used, and we reason that the conclusions of 
our investigation are not sensitive to this idealized treatment of the pumping. The constant 
that multiplies the density for the pumping is chosen separately for each cube so that 1 = 1 
for the strongest ray of the 128^ rays that emerge parallel to the magnetic field from the 
128^ grid points on the face of the particular cube. Representative computations also are 
performed when this pumping constant is chosen so that the strongest ray has 7 = 0.1 and 
7 = 10. The differences are not significant. 



4. Results 

When the effects of saturation are ignored (i.e., n is set to be AA/F regardless of the 
intensity) in solving equation (13) for the intensities, an unsaturated optical depth r can be 
obtained for the intensity 7, 

T = ln(7/7,) (24) 

Equation (13) is integrated along the paths of each of the 128^ rays that emerge from the 
grid points on the faces of the computational cubes. The unsaturated optical depth r at 
the peak intensity is then obtained for each ray. Representative histograms of these optical 
depths are given in Figure 1 for rays that are parallel and perpendicular to the direction of 
the average magnetic field in the cubes (with the normalization that 7 = 1 for the brightest 
ray in the parallel direction for each cube). The expected result that the parallel rays should 
have larger optical depths is evident for va/cs > 3. A similar trend with va/cs is evident 
in all of the 81 cubes that were studied. The 81 cubes were obtained by considering the 
solutions at nine equally spaced times (separated by At = 0.2L/cs) for each of the nine 
sets of initial conditions for which the MHD equations are integrated. Although there is 
considerable overlap in the histograms for va/cs = 3, the masing spots are represented by 
the brightest rays from the face of the cube and these are due to only a small number of the 
largest optical depths. There is no overlap even at va/ Cs = 3 for the tails of the distributions 
at large optical depth. It might be surprising that the perpendicular optical depths tend 
to be larger than the parallel optical depths in our calculations for Va/cs = 1. This result 
is somewhat artificial. It occurs because the turbulent velocities are smaller in the MHD 
computations for va/cs = 1 than for va/cs = 3 and va/cs = 10, and they also are smaller 
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than the thermal velocity Vth of the OH. When Vturb is less than vth in f± and /o, the effect 
of the MHD anisotropy of Vturh to enhance the parallel, relative to the perpendicular, rays is 
reduced or disappears. In contrast, the variation in (fiB/vth) along a ray is still significant (at 
least for the adopted ^Bavg/vth = 9) and reduces the optical depths for the a components, 
but not for the tt components since /o does not depend on B. We have verified that the 
histograms for the parallel and perpendicular rays coincide for va/cs = 1 when the Zeeman 
splitting is set to zero, as would be expected. 

A factor in the decrease in the anisotropy of the optical depths in Figure 1 with decreas- 
ing va/cs must be the decrease in the turbulent velocities that results when the initial va/cs 
is reduced while the initial rms turbulent velocities remain fixed in the MHD simulations. 
The rms values for the components of the turbulent velocities in a specific direction (i.e., 
parallel or perpendicular to the direction of the magnetic field) are approximately 0.9cs, O.lcs 
and OAcs for the turbulent cubes with initial va/cs —10, 3, and 1, respectively, in Figures 
1-4. However, this cannot be the only factor. In the turbulent cubes from the time evolution 
of the MHD equations that we use, the component rms velocities for va/cs = 3 decrease to 
approximately OAcg without a significant change in the anisotropy as indicated by histograms 
such as in Figure 1. A key difference emerges when we examine the dispersions (i.e., standard 
deviations) in velocity along the paths of individual rays. For the simulated turbulent cubes 
with va/cs — 3 and 10, the average dispersion in the component of the velocity along the 
line of sight is much smaller for rays that propagate parallel to the magnetic field than for 
rays that propagate perpendicular to the magnetic field. The two dispersions are equal for 
turbulent cubes computed with va/cs — 1 (for all va/cs, the rms of the components of the 
velocities and of the turbulent magnetic fields parallel and perpendicular to the magnetic 
field are equal). This indicates that the turbulent medium becomes significantly anisotropic 
when P = cl/v\ is near 0.1, at least in our calculations. 

The coupled equations (19) are now integrated with the effects of saturation included 
with rapid cross/velocity relaxation as described in the foregoing Section (i.e., n is given 
by equation 23) for turbulent cubes that are computed with va/cs = 3. Computations 
are performed for a range of intermediate angles as well as for propagation parallel and 
perpendicular to the magnetic field. At each angle, the coupled equations are solved for the 
rays that emerge from the 128^ grid points on the face of a cube. The peak intensity of 
the brightest ray at each angle is plotted in Figure 2 for representative results. To agree 
with the observations, all tt components and rays with large fractional linear polarization 
must be too weak to be detected. Hence, the focus is on the brightest rays. Since there is 
considerable scatter in the intensities in Figure 2, we encompass the intensity points with a 
shaded band to provide a sense of the behavior of the angular variation. The main result 
in Figure 2 is that the intensities within approximately 20° of the direction of the average 
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magnetic field are greater by a factor of about 100 or more than those at angles near 90°. 
All of these more intense rays near 0° are a components as indicated their high fractional 
circular polarization, as well as by their Doppler velocities which are quite close to -^[iBavg- 
All of the weaker rays near 90° are 100 % linearly polarized and are tt components (they are 
at essentially zero Doppler velocity). Note that, as discussed in §3, the pumping is chosen 
so that / = 1 for the brightest ray that is parallel to the average magnetic field in each 
cube. When there is no turbulence (the cube is uniform) and the masing is treated in the 
unsaturated limit, the angular variations of the intensities for a and tt components are given 
by the simple functions of equations (10) and (11). For comparison, these functions also are 
plotted in Figure 2 with tq chosen so that / = 1 at 0°, as well ( it also follows in this case that 
/ = 1 at 90°). Finally, not only is the polarization of the single brightest ray of a turbulent 
cube the same as that of a "pure" a component at angles less than about 20° (essentially 
100% circular), but the same is true for all of the rays at these angles that have intensities 
within 1% of the brightest rays. Likewise, the rays near 90° having intensities within 1% of 
the brightest rays at that angle are 100% linearly polarized (as expected for tt components) 
and have the Doppler velocity expected for tt components. Thus, only a components within 
about 20° of the average magnetic field will be detected from the turbulent cube in Figure 
2 when the observational sensitivity is 1% of the intensity of the brightest maser spot. 

A grayscale map of the intensities that emerge from the surface of a representative 
cube with vaIcs — 3 and parallel to the average magnetic field is shown in Figure 3. The 
general appearance seems compatible with that of observational maps. The spectra in the 
side panels of the three adjacent rays within a cluster demonstrates that a single feature of 
right circular polarization, a single feature of left circular polarization and a Zeeman pair 
can be in close proximity. The spectrum of the strongest feature in the entire map is a 
single feature that is 100% circularly polarized. All are at Doppler velocities expected for a 
components ( ±iJ,Bavg)- 

Finally, a histogram is given in Figure 4 for the ratio of the peak intensities of the weaker 
to the stronger feature of a Zeeman pair for all rays with intensities greater than 3% of the 
most intense ray for the two turbulent cubes with v^/cs = 3, one of which is the cube used 
in Figure 1. The 3% is an estimate for the weakest features that will be recorded in the 
observations. One of the histograms in Figure 4 demonstrates that the description of the 
turbulence being used here is consistent with the observational result that only one of the two 
circularly polarized Zeeman components frequently is detected for the mainline 1665/1667 
MHz masers. For the variations of the velocities and magnetic fields in other turbulent 
cubes (e.g., the other histogram in Figure 4), we find the full range of possibilities — in some 
the Zeeman pairs and unpaired components occur with comparable frequency, in others 
the Zeeman pairs are much more numerous and finally, there are many in which unpaired 
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components are completely absent. The rays in the histogram are subdivided according 
to peak intensity to explore whether the likelihood for an unpaired component depends on 
intensity. No such dependence stands out, though the Zeeman pairs are not among the most 
intense subgroup in this histogram. 

The most intense rays from a cube depend on the tails of the distributions in Figure 1 

at large optical optical depths. The tails of the distributions are, as is usually the case for 
tails of distributions, likely to be sensitive to details that are not well described in idealized 
models. However, the differences in the average optical depths in the distributions in Figure 
1 provide compelling and robust evidence that the a components will tend to be much more 
intense than tt components. 

5. Discussion 

For the MHD anisotropy to be effective in altering the optical depths, the turbulent 
velocities must be significant in comparison with the thermal dispersion of the OH molecular 
velocities. This minimum requirement is easily consistent with the observational data. The 
thermal velocity Vth for OH is ~ 0.3 km s^^. In contrast, velocity dispersions of at least a few 
km s~^ are typical within the region containing a cluster of maser spots (e.g.. Argon, Reid, & 
Menten 2000). The speed of sound for a gas of mostly molecular hydrogen at 100 K is 0.6 km 
s~^. A second condition — that the magnetic pressure exceed the thermal pressure (/3 < 1) 
— also seems to be satisfied. The observed magnetic field strengths typically are a few to 10 
milligauss. Hence, /? = c^/v\ ~ 0.1[(iV/10^cm"3)(T/100 K)/(i?/4mG)^] where a gas density 
of 10^ cm^^ H2 molecules is adopted as benchmark density based on calculations for the 
pumping of the masers (Cesaroni & Walmsley 1991; Pavlakis & Kylafis 1996). Although 
there are uncertainties in the number density and the kinetic temperature T, as well 
as in the precise maximum value of (3 that is allowed, it is reasonable to expect from the 
observations of B that va/cs > 3 for which our calculations in Figure 1 indicate that the 
magnetic field will be strong enough to create the necessary anisotropy. More generally, 
equipartition of energy between the turbulent magnetic field and the turbulent motions is 
suggested both by observations of the interstellar medium and by MHD simulations. Then, 
even in the absence of a significant uniform magnetic field, the turbulent magnetic field 
satisfies B/ (47rp)^/^ ~ Vt, which is equivalent to the statement /3 < 1 and that Va. > Q if the 
turbulent motions are supersonic. 

In addition to (1) the ratio of the turbulent to the OH molecular velocity and (2) the ratio 
of the thermal to the magnetic pressure (/?), another key consideration is the largest scale 
length (or minimum wavenumber) at which the turbulence is introduced into the medium 
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— the outer scale length. Since this determines the maximum length for an elongated 
masing volume, injection of the turbulence at the largest scale lengths is most effective for 
creating anisotropy in the optical depths. As noted in §3, a minimum wavenumber /Cmin 
that corresponds to the maximum wavelength L/2 is used for the calculations here. The 
relationship between the length L of the edge of the cubic volume within which the fields 
are computed and actual distances within interstellar clouds can be estimated from the size 
of a masing spot. In the maps that we compute such as Figure 3, the spectra and intensities 
change significantly over distances of only a few grid points. We naturally associate this 
dimension with the "spot size" ~ 10^^ cm for mainhne OH masers so that the 128 grid points 
along the edge of a cube correspond to L ~ 10^^ cm. This L is closer to the dimensions of 
one of the dozen or so clusters of masers in a masing cloud (e.g., Garcia-Barreto et al. 
1988) than to the entire cloud. A few computations with Vx/cg = 10 were performed where 
the maximum wavelength for injecting the turbulence was reduced to L/5. The differences 
between the mean optical depths parallel and perpendicular to the magnetic field analogous 
to those in Figure 1 are reduced, though the parallel optical depths can still be the larger by 
about 50% — sufficient for the maser emission along the magnetic field to be more intense 
by a factor of 100 or more than that perpendicular to B. 

The premise of our calculatons that the OH 1665/1667 GHz masers are not highly 
saturated is supported by the observation that the spectral line breadths tend to be narrower 
than the thermal breadths for the kinetic temperatures that are expected in the masing gas. 
The masers may still be somewhat saturated {R > T) as long as the rate 7 for velocity /cross 
relaxation satisfies 7 > i?. 

The OH masers in interstellar clouds are believed to arise in gas that has been infiuenced 
by shock waves. Hence, our use of decaying turbulence (as opposed to continuously driven 
turbulence) seems most appropriate. 

Dissipative processes will tend to erase structure in the gas at small scales. We need to 
check that the dissipation scale is smaller than our numerical resolution in the simulations 
and smaller than the maser spots ~ 10^^ cm. For N ~ lO^cm"^, T ~ 10^ K, and a fractional 
ionization x ~ 10~^, the dominant dissipative process will be ambipolar diffusion. In linear 
theory, ambipolar diffusion critically damps an Alfven wave when the wave frequency kvjy 
(calculated in the absence of ambipolar diffusion) is comparable with the collision frequency 
i^ni ~ xN{av)j^ between a neutral and an ion, where {cyv)j^^ = 1.9 x 10~^cm^s~^ (Kulsrud 
& Pearce 1969). This corresponds to Xd ~ 10^%T/100Ky/'^(3-^/'^{W cm-^/N){10-'^/x), or 
somewhat less than the size of a maser spot. 

The MHD anisotropy for optical depths parallel and perpendicular to the magnetic field 
also applies to other spectral lines — both masing and non-masing (or "thermal") spectral 



-18- 



lines. These differences in the optical depths may lead to anisotropic (or "m-dependent" 
) pumping as has been considered previously for the linear polarization of both masing 
(Western & Watson 1983, 1984) and thermal (Goldreich & Kylafis 1981, 1982; Deguchi & 
Watson 1984; Lis et al. 1988) spectral lines where the Zeeman splitting is much less than 
the spectral linebreadth. When applied to the pumping of OH masers, this anisotropy may 
also contribute significantly to the polarization of OH masers and to the suppression of the 
TT component. 

In contrast to our idealized calculation with phcnomcnological pumping rates and only 
two energy levels, Gray & Field (1995) calculated the excitation and transport of polarized 
OH maser radiation in a magnetic field including explicitly the magnetic substates of a 
number of the energy levels of the OH molecule. The important idea that a radiation 
would be strongly beamed in the direction parallel to the magnetic field and tt radiation 
beamed perpendicular to the magnetic field was recognized there. However, at that time 
no cause was recognized for suppressing the propagation of maser radiation perpendicular 
to the magnetic field in comparison with that parallel to the magnetic field, and hence for 
ultimately suppressing the tt components in comparison with the a components. 

Faraday rotation within a source tends to destroy the linear polarization of radiation and 

has sometimes been mentioned as a possible cause for the absence of the linearly polarized tt 
components and for the general tendency for circular polarization to dominate in OH masers. 
For unsaturated masing and large Faraday rotation, equations (10) and (11) would become 
(GKK, eq. 56 & 57) 

/± = (/e/2) exp[ro(l + |cos^|)V4] (25) 

and 

/o = (/e/2)exp[rosin2^/2] (26) 

The modification is thus equivalent to reducing the magnitude of the optical depth of the tt 
component by a factor of two at all angles while leaving the intensity of the a components 
essentially unchanged at angles near ^ = where they arc most intense. Its cfi^cct on the 
relative intensities of the tt and a components could be at least as great as what we have 
calculated for the influence of MHD turbulence. Faraday rotation depends on the column 
density of free electrons through the masing region and on the component of the magnetic 
field parallel to the fine of sight (e.g., Spitzer 1978, pg. 66). While an adequate estimate 
of the latter can be obtained from the observations, the column density of free electron is 
highly uncertain and reasonable estimates seem to allow the possibility that Faraday ro- 
tation could be cither large or negligible for the 1665/1667 MHz OH masers. However, a 
compelling argument against Faraday rotation as the general cause for the absence of the tt 
components is the observation that, although circular polarization dominates, linearly polar- 
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ized radiation frequently is detected in the a components. If Faraday rotation were always 
strong enough to suppress the tt components, it would also be strong enough to eliminate 
completely the linear polarization of the a components — in conflict with the observations. 
Again, considering Faraday rotation in the regime of saturated masing does not help. In the 
limit of large Faraday rotation and high saturation, the o components must be completely 
circularly polarized and, in principle, the three intensities /± and /q are equal at all angles 
(GKK, cq. 47). 

Although fractional linear polarizations of 10% or so are probably compatible with 
our calculations as they stand for suppressing the tt components, higher linear polariza- 
tions (> 30%) are detected in some clouds in a small fraction of the components and these 
are interpreted as a components (e.g., Garcia-Barreto et al. 1988; Gasiprong, Cohen, & 
Hutawarakorn 2002; Slysh et al. 2002). Note that circular polarization in a feature may not 
assure that the feature is a cr component since a fraction of the linear polarization can, in 
principle, be converted into circular polariztion by Faraday rotation or magnetorotation. If 
these are in fact a components, they present a challenge. This radiation must, presumably, 
be emitted at large angles Q relative to the magnetic field as seen from equation (13) and 
Figure 2. At these large angles, the optical depths of the tt components would be expected 
to be comparable to the optical depths of the a components and hence to be detected in 
similar (albeit small) numbers — in apparent conflict with the observations. However,we rea- 
son that there are additional considerations which may still allow the suppression of the 
TT components to be understood within the context of the anisotropy of MHD turbulence. 
(1) Our simpliflcations in considering only two energy levels and in ignoring any interaction 
between beams of maser radiation propagating in different directions may be inadequate. 
When these effects are included. Gray and Field (1995) find that the tt components are fur- 
ther suppressed if the optical depths favor directions along the field lines. (2) The anisotropy 
noted above in the optical depths of the radiation involved in the pumping will tend to cause 
the populations of the magnetic substates to be unequal, and also potentially contribute to 
suppressing the tt components. Despite this deficiency in the present calculations, it seems 
clear that MHD turbulence can create an anisotropic environment in interstellar clouds, and 
that the tt components will tend to be suppressed while the circular polarization is favored 
in the spectra of the mainline OH masers in this environment. 
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Fig. 1. — Representative histograms for the number of rays with unsaturated optical depths 
T. The optical depth is obtained at the peak intensity of each ray. The numbers are given 
as a fractions of the 128^ rays from the surface of a turbulent cube within intervals of 

0.5 in T. Separate histograms (as indicated) are shown for rays that propagate parallel 
and perpendicular to the average magnetic field of the turbulent cube for the three choices 
va/cs = 1, 3, and 10. 
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Fig. 2. — The intensity /(indicated by open circles and in units of the saturation intensity 
Is) for the brightest maser ray that emerges from a turbulent cube is plotted for several 
representative cubes with va/cs = 3 as a function of the angle between the ray and the 
magnetic field, and with the constant for the pumping chosen so that / = 1 at 0° for each 
cube. The intensities that individual a and tt components would have for unsaturated masing 
in a uniform medium are plotted (solid lines) for comparison (equations 10 and 11). At angles 
less than 55°, the fractional circular (V/I) and linear (Q/I) polarizations are plotted for this 
unsaturated a component (equations 12 and 13); for angles greater than 55°, the fractional 
linear polarization is plotted for the tt component (equation 15). The polarizations of the 
intensities of all rays at a particular angle that are computed using the turbulent cubes and 
for which the intensities are within 3% of the brightest ray at that angle have polarizations 
that are essentially identical with those of the unsaturated a and tt components in a uniform 
medium. The fractional circular and linear polarizations of the unsaturated components in 
a uniform medium are indicated by dashed and dotted lines, respectively. 
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Fig. 3. — Representative grayscale map of the intensities /(in units of Ig) that emerge from 
the face of a turbulent cube with va/cs — 3 and parallel to the direction of the magnetic 
field. The panels on the left show the spectra for the rays taken from selected locations 
on the map as indicated. The Doppler velocity scale is in thermal breadths Vth for the OH 
molecule. 
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Fig. 4. — Examples of the ratio of the weaker to the stronger of the a components of the 
Zeeman pairs for turbulent cubes with va/cs — 3. The lefthand histogram is computed 
for the turbulent cube in Figure 2. Only the rays with intensities greater than 0.03 of the 
brightest ray are included. The rays are grouped as fractions of the intensity of the brightest 
ray. Those with intensities from 0.03 to 0.1 of the brightest ray are indicated by black; those 
with intensities between 0.1 and 0.3 of the brightest, by gray; those with intensities between 
0.3 and the brightest, by white. The scale for each of these groups is normalized so that for 
each group, the sum in all columns of a type is one (and the sum of all columns is three). 
The bins are 0.05 in width. For some 75% of all rays in the lefthand histogram, the ratio is 
less than 0.05. The righthand histogram is an example where there are more Zeeman pairs. 



